Localized form of Fock terms in nuclear covariant density functional theory 
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In most of the successful versions of covariant density functional theory in nuclei, the Fock terms 
are not included explicitly, which leads to local functionals and forms the basis of their widespread 
applicability at present. However, it has serious consequences for the description of Gamow- Teller 
resonances (GTR) and spin-dipole resonances (SDR) which can only be cured by adding further 
phenomenological parameters. Relativistic Hartree-Fock models do not suffer from these problems. 
They can successfully describe the GTR and SDR as well as the isovector part of the Dirac effective 
mass without any additional parameters. However, they are non-local and require considerable 
numerical efforts. By the zero-range reduction and the Fierz transformation, a new method is 
proposed to take into account the Fock terms in local functionals, which retains the simplicity of 
conventional models and provides proper descriptions of the spin-isospin channels and the Dirac 
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Research on quantum mechanical many-body prob- 
lems is essential in many areas of modern physics. By 
reducing the many-body problem formulated in terms of 
A/'-body wave functions to the one-body level with local 
density distributions p(r), the density functional theory 
(DFT) of Kohn and Sham [l|] has accomplished great 
success. No other method achieves comparable accuracy 
at the same computational costs. Even though the ex- 
act form of the functionals is always extremely difficult 
to determine, this theory states that there exists a local 
single-particle potential l^s(^) so that the exact ground- 
state density of the interacting system can be reproduced 
by non-interacting particles moving in this local poten- 
tial. 

In nuclear physics, the DFT is used since the 1970s 
with great success to describe both ground states and 
excited states of nuclei throughout the periodic chart Q . 
In contrast to Coulombic systems, the spin and isospin 
degrees of freedom play an essential role, and the cor- 
responding functionals contain many phenomenological 
parameters whose numbers can be reduced considerably 
by the symmetries and properties of underlying theories. 
Therefore, the covariant density functionals [3|-[6| are of 
particular interest in nuclear physics. They exploit basic 
properties of QCD at low energies, in particular, sym- 
metries and the separation of scales They provide 
a consistent treatment of the spin degrees of freedom, 
they include the complicated interplay between the large 
Lorentz scalar and vector self-energies induced on the 
QCD level by the in medium changes of the scalar and 
vector quark condensates [8], and they include the nu- 
clear currents induced by the spatial parts of the vector 
self-energies. Furthermore the Lorentz symmetry puts 
stringent restrictions on the number of parameters in the 
corresponding functionals without reducing the quality 



of the agreement with experimental data. 

Over the years, a large variety of nuclear phenom- 
ena have been described successfully by these methods 
[El, B B [13 • the equation of state in symmetric and asym- 
metric nuclear matter, ground-state properties of finite 
spherical and deformed nuclei all over the nuclear chart, 
collective rotational and vibrational excitations, as well 
as fission landscapes and low-lying spectra of transitional 
nuclei involving quantum phase transitions in finite nu- 
clear systems. 

Most of these successes could only be achieved because 
the underlying functionals depend only on local densities 
and currents. The non-local Fock terms of the meson 
exchange potentials are not taken into account explic- 
itly and thus the theoretical framework is relatively sim- 
ple and computational costs are low. However, a few 
common problems have been found in the widely used 
relativistic Hartree (RH) method in the literature, 1) it 
is difficult to disentangle the effects of isovector-scalar 
and isovector- vector channels by the nuclear ground-state 
properties, unless a careful tuning is performed based 
on selected microscopic calculations [11|; 2) additional 
adjustment is needed for charge-exchange spin- flip ex- 
citations, such as Gamow- Teller resonances (GTR) and 
spin-dipole resonances (SDR); 3) nucleon-nucleon tensor 
interactions are practically missing so that the nuclear 
shell structures and their evolutions are not well repro- 
duced. 

Since a few years, relativistic Hartree-Fock (RHF) cal- 
culations are possible for spherical nuclei all over the nu- 
clear chart [12]. In contrast to the simple Hartree calcula- 
tions, they reproduce successfully the effective mass split- 
ting [12[, the spin-isospin resonances [13, 14], and shell 
structure evolutions [15|,[l6|. However, since RHF theory 
introduces non-local potentials VuFir^r')^ its theoretical 
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framework is much more involved and the computational 
costs are extremely heavy and the simplicity of the Kohn- 
Sham scheme is lost. This restricts the applicability of 
deformed RHF theory to a few very light nuclei [I^ and 
it prevents the treatment of effects beyond mean-field. 

It is therefore highly desirable to stay in nuclear physics 
too within the conventional Kohn-Sham scheme and to 
find a covariant density functional based on only local 
potentials, yet keeping the merits of the exchange terms. 
In this paper, we start from the important observation 
that one of the successful and most widely used RHF 



parameterizations PK02 [l6| includes only three rela- 
tively heavy mesons, a, cj, and p, but no pions. Since 
the masses of these mesons are heavy, the zero-range re- 
duction becomes a reasonable approximation. Moreover, 
with the help of Fierz transformation [l8|, |T9| , the Fock 
terms can be expressed as local Hartree terms. The main 
goal here is to investigate the validity of the local approx- 
imation to the exchange terms of the RHF theory. The 
isovector effective mass splitting and spin-isospin reso- 
nances will be focused, where the effect of the exchange 
terms is known to be significant [l2|-[l3 • 
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The basic ansatz of RHF theory is an effective Lagrangian density p^, |2| 
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in which nucleons are described by Dirac spinors that interact with each other by exchanging mesons and photons. 
Here, all the symbols have the same meanings as in Ref. [6[. With the Hartree- Fock and no-sea approximations, the 
system energy functional can be written as [2l| 



occ ^ occ ^ occ 

E={'i>o\H 1*0) =J2{a\a-p + pM\a) + -J2 {ab\ V{1, 2) \ba) - - ^ {ab\ V{1, 2) \ab) 



(2) 



ab 



ah 



where |<l>o) is the ground-state wave function as a Slater 
determinant, a, h represent the occupied single-particle 
states, and the two-body interaction ]/(!, 2) includes the 
following meson- nucleon interactions. 



K(l,2) = -[^,7o]i[^.7o]2l^a(l,2), 
K.(l,2) = [^c.7071i[^c.707m]21^c.(1,2), 
Fp(l,2) = [^p7o7^r]i-[^p7o7M^]2l^p(l,2), 
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and the photon- nucleon interactions. The corresponding 
finite-range Yukawa propagators of the mesons read 

A(n,r2) = - — ^ or = 2 , 2 W 
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with i = a^uj, p. 

At the scale of nuclear structure, the degrees of free- 
dom that have to be taken into account explicitly are 
pions and nucleons. The exchange of heavy mesons asso- 
ciated with short-distance dynamics cannot be resolved 
at low energies that characterize nuclear binding and can 
be represented by local four-point nucleon-nucleon inter- 
actions. This type of approaches has been demonstrated 
to be suitable for the descriptions of the nuclear ground- 
state and excited state properties [22|, [23| . In the heavy 
meson limit, ^ the Yukawa propagator in Eq. (|4]) 
can be expanded as 
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within the zero-order approximation. In next order, one 
finds the gradient terms, which are important for the 
description of ground-state properties like nuclear radii. 
They are not taken into account in the present calcu- 
lations of nuclear matter and collective excitations, be- 
cause the zero-order approximation is accurate enough 
for the physics discussed here. The zero-range coupling 
strengths a^^ related to the meson-nucleon coupling 
strengths and meson masses read 
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In present paper, the subscripts 5, T, PS^ and PV 
label the scalar, vector, tensor, pseudoscalar, and pseu- 
dovector channels, respectively, while without (with) t 
indicates the isoscalar (isovector) nature. 

According to the Fierz transformation [l8|, |l9| , the ex- 
change terms can be expressed as linear superpositions 
of direct terms with the equivalent coupling strengths 
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where j = {S,tS,V,tV,T,tT, PS.tPS, PV.tPV}, k 
{S^ V^tV}, and the transpose of matrix C reads 
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-8 -8 20 4 8 8 4 4 
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FIG. 1: (Color online) Zero-range coupling strengths of 
PK02-H in the isoscalar (left column) and isovector (right 
column) channels as a function of the baryonic density. DD- 
ME2 [2j parameters are shown for comparison. The dotted 
line is afs + afy 



In this way, the coupling strengths in ten channels are 
determined by the three density-dependent RHF parame- 
ters of PK02. The RHF equivalent parametrization thus 
obtained will be called PK02-H, and its zero-range cou- 
pling strengths are shown as a function of the baryonic 
density in Fig. [H They are compared with the corre- 
sponding strength parameters of the widely used density- 
dependent RH parametrization DD-ME2 [24], which in- 
cludes only V ^ and tV channels. In the following, 
we denote full Hartree-Fock results with parametriza- 
tion PK02 by RHF, the present results with PK02-H 
by RHF-loc, and the conventional Hartree results with 
parametrization DD-ME2 by RH. Very good agreement 
between both sets are found for the size and the density 
dependence of the two isoscalar channels S and V . This 
is not surprising because in the covariant framework the 
balance between these two channels is very delicate and 
plays an essential role in both mean-field and spin-orbit 
potentials. 

However, a substantial difference between RHF-loc and 
RH appears in the isovector channels. In DD-ME2 the 
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FIG. 2: (Color online) Dirac mass splitting between neutron 
and proton in neutron matter for RHF (shadow band) and 
RHF-loc (solid line) as the function of the baryonic density. 
The predictions of the Dirac- Br ueckner- Hartree-Fock calcula- 
tion [2^ are shown as close symbols. 



isovector-scalar ^-meson is not explicitly included and the 
isovector-vector p-meson determines the full isospin de- 
pendence. Eq. ([7|) shows that the W term in PK02-H is 
strongly repulsive. It is modified by the strongly attrac- 
tive iS term and for the large components we find as in 
the isoscalar case a delicate balance between scalar and 
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vector channels, 
(see the dotted line in Fig. [1]). 

The scalar fields in the isovector channel determine 
the Dirac mass splitting between neutron and proton in 
asymmetric nuclear matter. Due to the absence of the tS 
channel, it vanishes identically for RH. This contradicts 
to the results of microscopic Dirac-Brueckner-Hartree- 
Fock (DBHF) calculations [25|. In Fig.O the Dirac mass 
splittings in neutron matter are presented for RHF-loc 
(solid line) and DBHF [25] (close symbols). The corre- 
sponding cuves are compared also with results for the 
full RHF calculations. In RHF, the Dirac mass is mo- 
mentum dependent and therefore as a shadow band is 
shown, where its upper (lower) border corresponds to the 
results calculated with neutrons at zero (Fermi) momen- 
tum. Note that the proton Fermi momentum vanishes in 
neutron matter. In the present RHF equivalent calcula- 
tions, the Dirac mass does not depend on the momentum, 
but it maintains the characteristic behavior of the RHF 
results. It is also in very good agreement with the mi- 
croscopic DBHF results. This indicates the constraints 
introduced by the Fock terms of the RHF scheme into the 
tS channel of the present density functional are straight 
forward and quite robust. 

In panels (e)-(j) of Fig.[Tl it is found that the remaining 
six channels of PK02-H do not vanish. They are explic- 
itly determined in Eq. ([7]) by the exchange effects of the 
RHF scheme. In contrast, the coupling strengths ce^^^p^ 
and <^[^)py are practically out of control in usual fit- 
ting processes since they are absent in the description of 
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FIG. 3: (Color online) Upper panel: Strength distribu- 
tion of the GTR in ^°^Pb calculated with a Lorentzian 
smearing parameter of 1 MeV. The experimental peak en- 
ergy 26] is denoted by an arrow. Lower panel: Diagonal 
matrix elements of particle-hole residual interactions for the 
(z/l/iii/2)~"^(7rl/i9/2)^ configuration. Further details are given 
in the text. 

ground states due to parity conservation. However they 
play a crucial role in the properties of charge-exchange 
modes. 

In the following, we demonstrate their importance 
for Gamow- Teller and spin-dipole resonances excited by 
charge-exchange reactions in spherical nuclei. Based 
on the unperturbed single-particle spectra calculated 
by RHF theory, we perform calculations based on the 
proton-neutron random phase approximation (RPA) us- 
ing the particle-hole {ph) interactions 

^t5(l,2) = afs[jor]i ■ [jorh^iri - r2), 
l/ty(l,2) = afy[7o7^r]i • [70 2(^(^-1 - ^2), 
VtT{h2) = af^[7oa^"r]i • [7o^M-r]2(^(ri - r2), (8) 
1/^P5(1,2) = afpsilolbr]! • [7075^] 2 (^(r^i - ^^2), 
l/tpy(l,2) = afpy[jol5l^r]i • [7075 2(^(^-1 - ^^2), 

while all isoscalar contributions vanish. 

In the upper panel of Fig. [H the GTR strength distri- 
bution in ^^^Pb for the RHF equivalent RPA is compared 



with the results of RHF+RPA and RH+RPA For 
RH-FRPA, a phenomenological Migdal term had to be 
introduced in order to reproduce the GTR data in ^^^Pb 
p7j . Here the results with and without this Migdal term 
are shown. The present RPA calculations lead to quite 
similar results as those by the original RHF+RPA [l3| 
calculations. The slight difference is due to the zero- 
range approximation in Eq. (|5j). The black arrow indi- 
cates the experimental peak energy [26| . It is reproduced 
very well for RHF and RHF-loc, and it is fitted for RH 
by adjusting the Migdal term. 

To understand the physical mechanisms, we present in 
the lower panel of Fig. [3] the diagonal matrix elements of 
an important GT ph configuration {i^lhii/2)~^{'7^^hg/2)^ • 
The total matrix element (full) is decomposed in various 
ways: For the original RHF+RPA [13] calculations, the 
total ph strengths is essentially determined by the deli- 
cate balance between the a- and cj-mesons contributing 
via the exchange term to the isovector channel. The orig- 
inal p-meson gives only a minor contribution. The effects 
of a- and cj-mesons are repulsive and attractive, respec- 
tively, since they contribute only via the Fock terms. 

A very different decomposition is shown for RHF- 
loc, where the matrix elements in the five channels 
{tS^tV^tT^tPS^tPV) are compared. Of course, accord- 
ing to Eq. (O, all three mesons contribute here. It is 
seen that the major contribution comes from the re- 
pulsive tT channel. It is reduced by the tPV chan- 
nel. This can be understood by the following analysis. 
For the charge-exchange spin- flip modes, the dominant 
residual interactions are represented by the matrix el- 
ements of the operator [ar] • [ar] in the large compo- 
nents of Dirac spinors. From Eqs. (|8]), one can readily 
find that they are composed of [a'^^r] • [(Jijr] ex and 
[757*r] • [757ir] oc —afpy. This decides the signs of these 
contributions and their amplitudes roughly differing by 
a factor of 4. Furthermore, the net contribution is then 
approximately proportional to 

2afT-afpy = -^{a^ + a^), (9) 

according to the intrinsic relations in Eq. ([7]). In other 
words, the total ph strengths are determined by noth- 
ing but the delicate balance between the S and V chan- 
nels (corresponding to the a- and cj-mesons), which is 
well established in nuclear covariant density functionals. 
Therefore, the constraints introduced by the Fock terms 
of the RHF scheme into the ph residual interactions are 
also straight forward and quite robust. 

In the case of RH, the contributions from the p-meson, 
the Migdal term, and the pseudovector coupling of the 
pion are shown. It is found by fitting to the data that the 
free pion residual interaction is attractive and the Migdal 
term is repulsive. 

The local representation of the Fock terms can be fur- 
ther examined by another important charge-exchange 
spin- flip mode, the spin-dipole resonance (SDR). In 
Fig. m the strength distributions of SDR in ^^Zr of the 
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FIG. 4: (Color online) Strength distributions of the SDR in 
^°Zr calculated by (a) RH, (b) RHF-loc, and (c) RHF with a 
Lorentzian smearing parameter of 2 MeV. The dash-dotted, 
dotted, dashed lines show the = 0~, 1~,2~ contributions, 
respectively, while the solid lines show their sum. The exper- 
imental peak energy [2^ is denoted by the arrows. 



RHF equivalent RPA are compared with the results of 
RHF+RPA and RH+RPA with the Migdal term. The 
experimental peak energy [28| can be excellently re- 
produced by the original RHF+RPA and the present 
RHF equivalent RPA calculations. However, the tradi- 
tional RH+RPA calculations present a more fragmented 
structure due to the poorly adjusted isovector ph resid- 
ual interactions [29|. This is remarkably improved in 
the present RHF equivalent functional. Therefore, not 
only the total strengths but also the individual contri- 



bution from different spin-parity components are al- 
most identical to the corresponding RHF+RPA results. 
In particular, this also applies for the energy hierarchy 
E{2-) < E{1-) < E{0-). 

The above investigations on charge-exchange excita- 
tions illuminate that the success of the RHF+RPA ap- 
proach has a solid foundation. The essential roles played 
by the a- and cj-mesons via the Fock terms can be suc- 
cessfully folded into the equivalent interactions in the 
Hartree level. This achievement makes the theoretical 
framework much simpler and the computations much 
more economical. 

In summary, a new local RHF equivalent covariant 
density functional is proposed. By the zero-range reduc- 
tion and the Fierz transformation, the constraints intro- 
duced by the Fock terms of the RHF scheme are fully 
taken into account. In this way, the advantages of exist- 
ing RH functionals can be maintained, while the prob- 
lems in the isovector channel can be solved. This opens a 
new door for the development of nuclear local covariant 
density functionals with proper isoscalar and isovector 
properties in the future. As a consequence, the results 
of this investigation point in the direction of the con- 
ventional Kohn-Sham picture to represent the non-local 
exchange terms in the form of local potentials. 
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